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Abstract 

The strong K~p scattering length is extracted within chiral SU(3) unitary approaches 
from a very large variety of fits to low-energy K~p scattering data. Very good overall 
agreement with available scattering data is obtained and the resulting scattering length is 
compared with the new accurate kaonic hydrogen data from DEAR. The pole structures 
of the obtained fits to experiment are critically examined. 
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1 Introduction 



The approximate chiral symmetry of the Hght quarks -u, d, s in quantum chromo dynamics plays 
a crucial role in describing low-energy hadronic processes. Chiral SU(3) symmetry is exact in 
the limit of vanishing light quark masses, but finite quark masses, m„, 777.^,777.5, induce explicit 
symmetry- breaking. While for the very light quarks u, d the symmetry- breaking corrections 
are in general small, it remains unclear whether the strange quark is still light enough to be in 
the chiral regime. 

In this respect, the KN system provides a good testing ground for chiral SU(3) dynamics 
and the role of explicit chiral symmetry-breaking due to the strange quark mass. High-precision 
K~p threshold data, such as K''p scattering, the ttE mass spectrum and the precisely measured 
K~p threshold decay ratios set important constraints for theoretical approaches. Recently they 
have been supplemented by the new accurate results for the strong interaction shift and width 
of kaonic hydrogen from the DEAR experiment [1] which reduced both the mean values and 
error ranges of the previous KEK experiment [2] . 

The existence of the A (1405) resonance in the K~p channel just below its threshold makes 
the loop expansion of chiral perturbation theory (ChPT) inapphcable. In this regard, the 
combination of ChPT with non-perturbative coupled-channel techniques based on driving terms 
of the chiral SU(3) effective Lagrangian has proven useful by generating the A(1405) dynamically 
as an I = KN quasibound state and a resonance in the vrS channel. 

However, the recent and precise DEAR measurement appears to be in disagreement with 
the K~p scattering length derived from scattering data as pointed out in [3]. A thorough 
investigation of the low-energy K~p interaction within chiral unitary approaches has reinforced 
the question of consistency of the DEAR experiment with K~p scattering data [4]. In contrast, 
it was claimed very recently in [5, 6] that within a chiral unitary approach both the scattering 
and the DEAR data can be accommodated. 

It is evident from these recent investigations that the K~p system remains a topic of great 
interest and is under lively discussion, see also [7, 8] . The aim of the present work is to shed 
some more light on this issue by providing a conservative range for the K~p scattering length 
constrained solely from K~p scattering data. This is accomplished within different variants of 
chiral unitary approaches and by performing a very large number of fits to experiment, in order 
to reduce the inherent model-dependence of these frameworks. The obtained realistic range for 
the K~p scattering length is then compared with the one derived from the DEAR and KEK 
experiments. 

This work is organized as follows. In the following section we present the basic formalism. 
Our results for K'p scattering are shown in Section 3. Section 4 contains an overview of the 
pole structures of the obtained solutions in the A(1405) region. We summarize our findings 
in Sec. 5. Some of the detailed results of this investigation are collected and displayed in the 
appendix. 

2 Formalism 

In this section we illustrate the underlying features of the approach. Details of the formalism 
have already been presented in previous works [4] and will not be repeated here. The starting 
point is the chiral effective Lagrangian C — jCij) + jC^b which describes the coupling of the 



2 




(a) (b) (c) (d) 

Figure 1: Shown are the leading order contact interaction (a), the direct (b) and crossed (c) 
Born term as well as the next-to- leading order contact interaction (d). Solid and dashed lines 
represent baryons and pseudoscalar mesons, respectively. 



pseudoscalar meson octet {n,K,rj) to the ground state baryon octet {N, A,I1,E). Both the 
purely mesonic piece >C<^ and the meson-baryon Lagrangian jC^b are employed up to second 
chiral order. 

Due to the nearby A (1405) resonance unitarity effects from final state interactions are im- 
portant for KN scattering and must be included in a non-perturbative fashion. To this aim, the 
relativistic effective Lagrangian is utilized to compute the tree level amplitude T^6,ia(s, cr, a') 
of the meson-baryon scattering processes 0j-Ba — > (pjB^ (with spin indices a, a') at invariant 
energy squared s. This amplitude is the driving term in the coupled-channels integral equation 
determining the meson-baryon T-matrix. 

In the literature, the effective meson-baryon Lagrangian has been used as interaction kernel 
at different levels of sophistication. While only the Weinberg- Tomozawa term from the covariant 
derivative is taken, e.g., in [9], the Born terms are included in [10]. In [11], on the other hand, 
the Lagrangian of second chiral order is added which yields additional contact interactions. 
In order to provide an estimate of the model-dependence of such approaches, we will discuss 
the following three different choices for the amplitude Vjh,ia{s,fl;a,a'). First, only the leading 
order contact (Weinberg-Tomozawa) term is taken into account, see Figure la. Subsequently, 
the Born diagrams are included, see Figs, lb and c. In the third approach we add the contact 
interactions from the Lagrangian of second chiral order, (Fig. Id). For brevity, we will 
refer to these variants as "WT" (Weinberg-Tomozawa), "WTB" (Weinberg-Tomozawa -|- Born 
diagrams) and "full" (including also the higher order contact terms), respectively. 

For each partial wave I unitarity imposes a restriction on the (inverse) T-matrix above the 
pertinent thresholds 

ImT.-i = (1) 

with c[cm being the three-momentum in the ccntcr-of-mass frame of the channel under consid- 
eration. Since we are primarily concerned with a narrow center-of-mass energy region around 
the KN threshold, it is sufficient to restrict ourselves to the s-wave (matrix) amplitude V{s). 
We write the inverse of the T-matrix as (suppressing the subscript I {— 0) for brevity) 

T-i = y-i + G , (2) 

which yields after inversion 

T =[1 + V ■G]-^V. (3) 
The quantity G is the finite part of the scalar loop integral G 

^^^'^ " J {2ttY [{q - p)2 - M| + te] \p^-ml + le] ' 
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where Mb and are the physical masses of the baryon and the meson, respectively. In 
dimensional regularization one obtains for the finite part G of G [4, 10], 



G{q') = a(/x) + 
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(m^ + Ms)2 - q 




(5) 



where is the regularization scale. The subtraction constant a(/x) cancels the scale depen- 
dence of the chiral logarithms and simulates higher order contributions with the value of a(//) 
depending on the respective channel. Stated differently, this introduces some additional SU(3) 
breaking beyond the use of the physical masses in the kinematics and loop functions. Eq. (3) 
is a matrix equation with the diagonal matrix G collecting the loop integrals in each channel. 
This amounts to a summation of a bubble chain to all orders in the s-channel, equivalent to 
solving a Bethe-Salpeter equation with V as driving term, where all momenta in V are set to 
their on-shell values. This so-called on-shell scheme reduces the full Bethe-Salpeter equation to 
the simple matrix equation (3). (Note, however, that in the presence of the crossed Born term 
(Fig. Id) this simplification must be treated with care due to the appearance of unphysical 
subthreshold cuts [4].) 

Moreover, the Coulomb interaction has been shown to yield sizable contributions to the 
elastic K~p scattering amplitude up to kaon laboratory momenta of 100-150 MeV/c [14]. We 
have thus taken into account the electromagnetic interactions as well, cf. [4] for details. 



3 Results 

In this section, we present our results for K~p scattering and the resulting prediction for the 
K~p scattering length, ax-p- Low-energy antikaon-nucleon scattering and reactions have been 
studied experimentally decades ago [15-20]. The available data (admittedly with large errors) 

are mostly restricted to K~ momenta above 100 McV/c. Further tight constraints are imposed 
by the accurately determined threshold branching ratios into the inelastic channels ttS and 
7r°A [21,22]: 

,.£||^-l4|l. 2.3a ±0.04. 

1 (A p — > TT L + ) 

T[K p — > all inelastic channels) 

Rn = ^^J(^-P--'^^ ^ ^ 0.189 ± 0.015 , (6) 
r(i\ p neutral states) 

and by the ttE invariant mass spectrum in the isospin / = channel [23]. 

In our analysis we have restricted ourselves to pure meson-baryon scattering and have not 
included the processes K~p — > 7A(1405) [24], 7p — > K*A{1A05) [25] which is now experimen- 
tally under investigation at SpringS/Osaka [26] and, in the near future, also at ELSA (Bonn), 
and K^p 7r°7r°S° [27] which is already measured [28]. Reactions including the coupling to an 
external photon such as K~p — > 7A(1405) and jp —>■ iir*A(1405) require substantial extension 
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of the chiral unitary approach apphed here as illustrated, e.g., in [29], whereas the three-body 
final state in K~p — > 7r°7r°E° introduces additional model-dependence. Hence, these processes 
are beyond the scope of the present investigation. 

Our approaches have six subtraction constants a(/i) in the different channels and up to 
eight parameters which are varied in the fits within generous limits: the decay constant / and 
in the full approach the higher order couplings 6j, dj. The axial vector couplings D, F which 
enter the Born diagrams (Fig. lb and c) are kept fixed at the values D — 0.80, F — 0.46 
extracted from semileptonic hyperon decays [30]. We have purposely chosen ample ranges for 
the parameters of our approach, in order to be able to take into account a large variety of 
qualitatively different fits to K^p scattering data. Since our concern here is to predict the 
strong K^p scattering length only from K~p scattering data, we do not impose additional 
phenomenological constraints, e.g., the analysis of [31] which includes r] photoproduction as a 
high quality data set. In addition, the framework chosen in [31] does not exactly coincide with 
any of the approaches considered in the present work. We can therefore not expect the same 
values for the coupling constants. 

We perform an overall least-squares fit to available low-energy K~p scattering data for 
the three different approaches, "WT", "WTB" and "full". To this end, we first calculate the 
individual for the i-th observable and divide by the number of pertinent data points Ui. The 
total per degree of freedom (d.o.f.) is then defined as [32] 

^ Ei^» X| 
d.o.f. N{Y.ini-p)^ nr 

where N is the number of obscrvablcs and p the number of free parameters in the approach. 
This definition of the x^/d.o.f. function generalizes the standard x^/d.o.f. for a single observable 
and has the advantage that all observables are weighted equally regardless of the number of 
data points. If one were to use instead the definition x^/d.o.f. = EjX?/(Ei'^j ^ p) iii which 
all data points from different observables have the same weight then single-valued observables 
(such as branching ratios) would be dominated by observables with many data points (such as 
scattering data). Note that the definition in Eq. (7) reduces to the latter one if all observables 
have the same number of data points. 

With respect to the work [4] we have significantly improved our fitting procedure which now 
involves the combination of a Monte Carlo routine with a conjugate gradient method [33]. This 
allows us to perform a large number of different fits to data distributed in parameter space so 
that the model dependence of the results is reduced and a realistic error range for the K^p 
scattering length derived from scattering experiments can be provided. 

We obtain a large number of fits which describe the low-energy K^p scattering data very 
well with the minimum x^/d.o.f. values given by 1.28, 0.88 and 0.71 for the different approaches 
"WT", "WTB" and "full", respectively. In general, if is a function of n parameters 6 = 
[6i, . . . , On) the standard error range of these parameters is given by the condition [34] 

x'(0)=xLn + Ax^ (8) 

and Ax^ is derived from the p-value of the x^ probability distribution function with the per- 
tinent number of degrees of freedom. Strictly speaking, this relation only holds if the method 
of least squares is applied to one single experiment and the associated fit function depends lin- 
early on the parameters d. In the present investigation the situation is more involved: the free 
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parameters of the approach enter in a highly non-hnear way into the calculation of observables 
and the x^-function which is minimized combines a variety of measurements of different quan- 
tities. Nevertheless, we will adopt the standard definition of a confidence region, Eq. (8), since 
one can expect it to be a reasonable approximation — at least in the vicinity of the minimum 
of the x^-function, where its shape should be nearly parabolic. Our fit includes a total number 
of 171 data points and either 7 ("WT", "WTB") or 14 ("full") parameters. The region of one 
standard deviation (i.e. 68.27% confidence level) is then found by adding Ax^/d.o.f. — 1.05 to 
the minimal x^/d.o.f. as defined in Eq. (7) regardless of the approach ("WT". "WTB" or "full") 
since the difference in the number of parameters causes only tiny modifications in A^^/d.o.f.. 

As it turns out, we also obtain fits that although having a relatively low overall x^/d.o.f. 
fail miserably in one or two observables and can thus not be classified as "good fits". In 
addition to the bound on the overall x^/d.o.f. we thus demand that each individual observable 
be reproduced with a xl/i^i value which does not exceed the overall x^/d.o.f. by more than a 
factor of four. Fits which do not meet this additional criterion are grouped instead with fits with 
the lowest x^/d.o.f. value which satisfies this constraint. This specific choice has proved useful 
in practice. These two goodness-of-fit criteria, a bound on the overall ^^'^ on the individual 
Xi, determine the error regions specified in the following for all parameters and observables, 
while the presented central values correspond to the fits with minimal x^- We also investigate 
analytic properties of the fits, in order to sort out solutions with unphysical pole structures. 
This issue will be discussed in detail in the next section. 

The scattering length ax-p is given by the strong interaction T matrix at threshold 

dK-p = Tx-p-^K-pis) 

The values corresponding to the best fits in the three approaches are 

"WT" : aK-p = (-0.73 + i 0.91) fm , 

"WTB" : ttK-p = (-1.09 + i 0.84) fm , (10) 
"full" : aK-p = (-1.05 + i 0.75) fm . 

The errors of the real and imaginary parts are, of course, correlated and the la regions in the 
complex ttK-p plane are depicted in Fig. 2. While the absolute minimum of x^/d.o.f. is the 
lowest in the full approach, the x^ function rises steeper in the "WT" and "WTB" approaches 
leading to a smaller la confidence region in the ax-p plane. 

We also extract the KN s-wave scattering lengths a^, in the isospin limit of equal up- and 
down-quark masses. Since we neglect isospin-breaking corrections in the Lagrangian from which 
the interaction kernel of the coupled-channels calculation is derived, taking the isospin limit 
amounts to replacing the physical masses of the particles that enter all kinematic quantities 
and the scalar loop integrals G by a common mass for each isospin multiplet and we disregard 
any electromagnetic effect. For the different approaches we obtain the central values 

"WT": ao = (-1.45 + i 0.85) fm , ai = (0.65 + i 0.76) fm , 

"WTB": ao = (-1.72 + i 0.77) fm , ai = (0.09 + i 0.76) fm , (11) 

"full": ao = (-1.64 + i 0.75) fm , ai = (-0.06 + i 0.57) fm . 

In Fig. 3 we show the error ranges and compare the results with values found in similar chiral 
unitary approaches [5,6,9,10] and a multichannel dispersion relation analysis of KN scattering 



.s=rm,„_-l-M„~|2 



(9) 



6 



"WT" 



WTB^ 



1.2 
1.0 



t< 0.8 



0.6- 



0.4 




1.2 

1.0 



^ 0.8- 



0.6- 




0.4 



-1.5 -1.2 



-0.9 -0.6 -0.3 0.0 

Reax-p (fm) 



0.3 



-1.5 -1.2 



-0.9 -0.6 -0.3 0.0 
Reax-p (fm) 



0.3 



"full" 




-1.5 -1.2 -0.9 -0.6 -0.3 0.0 0.3 

Reax-p (fm) 

Figure 2: Real and imaginary parts of the scattering length ax-p for the three approaches. The 
circles indicate the result of the best fits, the shaded areas represent the la confidence region. 



7 



"WT" 



'WTB^ 



1.2 



0.8 



a 0.6 



0.2 



' 1 ' 1 
V 


1 1 1 1 1 1 1 


n / \ 




' ( ) 

: Loop/ 








zv-^ — ^ 






O *\ 7 


♦ ■ 


□ 




♦ 




1 


-2.0 -1.5 -1.0 -0 


.5 0.0 0.5 1.0 1.5 



ReOj (fm) 



1.2 



0.6 



0.2- 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 



' 1 ' 1 
V 


1 1 1 1 1 1 1 1 


- 

- / clqi 


^-..^^^ 


- / 




- \ A / 


J 




Vv * 


♦ ■ 


□ 




♦ 





1 



-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

ReOj (fm) 



"full" 



_ 1 ' 1 ' 1 


' 1 ' 1 ' 1 ' 1 


□/ \ 




- / flo J 




- \ A / 








♦ ■ 






♦ 




1 


-2.0 -1.5 -1.0 -C 


.5 0.0 0.5 1.0 1.5 



ReOj (fm) 

Figure 3: Real and imaginary parts of the isospin KN scattering lengths Oq and ai for the 
three approaches. The circles indicate the result of the best fits, the shaded areas represent the 
la confidence regions. Our results are compared to the values found in [35] (triangle up), [9] 
(star, the value ao = (—2.24 + 1.94i)fm is not shown), [10] (triangle down), [5] (filled diamond 
for fit y4j, empty diamond for fit B^), [6] (filled square for fit I, empty square for fit 11). 



data [35]. One observes consistency of all three approaches ("WT", "WTB", "full") within error 
bars and agreement with most of the values from previous investigations. However, our results 
do not agree with fit in [5] (as already discussed in [7]) and the similar fit I in [6]. They also 
disagree with the ai value in [9], where a variant of the "WT" approach was employed. The 
approach utilized in [10] nearly coincides with "WTB" in the present work, however only one 
common subtraction constant for all channels was employed in [10], while in the present work 
we have the freedom to vary six (isospin symmetric) subtraction constants. This explains why 
the imaginary part of oq in [10] is larger and outside the la range of the present calculation. 

Due to the very high statistics of the present investigation our results provide a realistic 
error range for the isospin scattering lengths oq and Oi within chiral unitary approaches. As 
shown in [36] these quantities are an important input in the theoretical analysis of the upcoming 
spectroscopy study of kaonic deuterium [37] and further anticipated experiments with even more 
complex light kaonic nuclei [38] at DA$NE. 

In the presence of electromagnetic corrections the ground state strong energy shift AE" and 
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width r of kaonic hydrogen are related to the K p scattering length ax-p via [3] 

AE-^r^ -2a^pLlaK-p [1 - 2Q;//,(ln a - l)aK-p] , (12) 

where Hc is the reduced mass of the K~p system and a is the fine-structure constant. The 
obtained predictions for AE and F, which are solely based on K~p scattering data and the vrS 
invariant mass spectrum, are presented in Fig. 4 for the three different approaches. They are 
compared to the recent experimental determination of kaonic hydrogen observables by both the 
KEK [2] and the DEAR [1] collaborations. Figure 4 constitutes together with Figs. 2 and 3 the 
main result of the present work. The shaded areas in the plots of Fig. 4 represent smoothened 
areas which correspond to different upper limits of the overall x^/d-o.f. (with the additional 
constraint that each individual observable is reproduced by the fit with x^/ui of at most four 
times the upper limit of x^/d.o.f.) and are drawn on the basis of several thousands of fits (e.g. 
more than 7000 for the full approach) . We point out that a comparable statistical exploration 
of parameter space in chiral unitary approaches for interactions has not been attempted 
before; it provides for the first time a realistic estimate of theoretical uncertainties within this 
framework. 

Regardless of the chosen approach the fits with minimal overall agree nicely with the 
result of the KEK experiment, while the la confidence region, which is bordered by the dashed 
line in the plots of Fig. 4, has no overlap with the error ranges given by the DEAR experiment. 
As discussed above the standard definition of the la confidence region by means of Eq. (8) is 
not strictly applicable in the present investigation, where the fit function is non-linear in the 
parameters and the fit incorporates a variety of different observables. Therefore we refrain from 
showing 2a and 3(j confidence regions since application of the standard error estimation seems 
more questionable in these cases. 

Instead we plot regions that correspond to quadratically increasing upper limits of the 
overall x^. Note that fits which are compatible with the error ranges given by DEAR have an 
overall xVd.o.f. of at least 6.1, 5.5, 3.3 in the "WT", "WTB", "full" approach, respectively, 
with elastic K~p scattering being the largest source of disagreement. Finally, for the best fit of 
each approach we present in Table 1 the numerical values of the fitted low energy constants /, 
6o, i>£), bp, di, d2, d^, d^ and the subtraction constants in the loop integrals G. Note that in the 
fits we have allowed for broad ranges for the subtraction constants a^B- In fact, in the "WT" 
and "WTB" approaches the resulting axs are roughly one order of magnitude larger than the 
remaining subtraction constants. However, the fits are not very sensitive to variations in this 
parameter such that the value is only slightly increased if axs is reduced to the same size 
as the other a^^. The detailed comparison of experimental input from KN and ttE scattering 
with our fits is compiled in the appendix. 

4 Resonance poles 

In order to cover a substantial region in parameter space we put only very loose constraints 
on the numerical values of the low energy parameters of the chiral effective Lagrangian, i.e., 
we solely fix the order of magnitude as motivated by the naturalness assumption of couplings 
in the effective field theory. The subtraction constants in the loop integrals G are permitted 
to vary in even larger ranges. Starting from randomized initial values of the parameters the 
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Figure 4: Strong energy shift AE and width F of kaonic hydrogen for the three approaches. 
The shaded areas represent different upper limits of the overall x^/d.o.f. The la confidence 
region is bordered by the dashed hue. See text for further details. 
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Table 1: Numerical values of the fitted couplings and subtraction constants corresponding to 
the best fits in the three approaches. The empirical value of the average meson decay constant 
is / ~ 100 MeV. 
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Monte Carlo routine utilized in the present investigation generates a vast number of fits and for 
all of them the poles of the T-matrix in the complex W — plane are determined. Although 
we work in the physical basis, where isospin is broken by the physical masses of the particles 
in the loop integrals, we can classify the poles as being mainly of isospin / = 0, 1,2 by their 
impact on channels (or channel combinations) which contain only one isospin component (e.g. 
rjA, r/S^). The pole positions of the fits will serve as an additional constraint to rule out certain 
fits which must be considered as unphysical as we will explain in the following. 

While resonances are generally associated with poles in unphysical Riemann sheets of the 
complex W plane, the solutions should be free of poles in the physical Riemann sheet as required 
by the postulate of maximal analyticity. Since it is not possible to numerically explore the entire 
upper half-plane, we must choose a finite region to search for poles in the physical sheet. To 
this end, we dismiss all fits which exhibit a pole at a distance of less than 250 MeV from the real 
axis in the relevant energy region (1.25 - 1.50 GeV). This selection criterion ensures that, even 
if such pathological poles were to exist, their influence on the real axis would be negligible. In 
particular, it guarantees that the Wigner bound, which is based on causality and sets a lower 
limit for the derivative of the phase shift with respect to energy, is not violated [7,39]. Note that 
in a similar manner the Wigner condition has been employed in [40] to constrain the numerical 
values of parameters in the context of unitarized chiral effective field theory. 

Secondly, we reject fits which have a resonance pole on the relevant unphysical Riemann 
sheet that is located less than 2.5 MeV below the real axis and thus corresponds to a resonance 
with a width of less than 5 MeV. Such resonances would be one or two orders of magnitude 
narrower than what one would expect from the characteristic time-scale of the strong interac- 
tions of about 10~^^ s. Lifetimes of this order correspond to widths of several tens to hundreds 
of MeV, in agreement with typically observed hadronic resonances in this energy region. At 
present, there is no experimental indication for an exotically long-lived stated in the energy 
interval under consideration and we can safely ignore such fits. 

The third item concerns the ttE event distribution [23] which clearly shows a peak corre- 
sponding to the A(1405) resonance. We adopt the approach advocated in [10] which describes 
the experimental vr'E^ event distribution as originating from a generic / = source made up 
of unknown shares of vrS and KN states. As it happens we observe fits which do not exhibit a 
true resonance structure around 1.4 GeV, but merely show a broad bump that is generated by 
the intricate superposition of the two source states, ttE and KN. As the ttE event distribution 
is not normalized, the normalization constant in the fit can be tuned in such a way that the 
of the pertinent fit can have a relatively low value. However, these fits do not have an isospin 
zero pole on the unphysical sheet at a position which could be associated with the A(1405). In 
fact, if the ttS event distribution were simply approximated by the invariant mass distribution 
of / = ttE states, see e.g. [9,41], one would observe no peak structure at all for these fits. 
Taking the well-estabhshed four star resonance A(1405) for granted, one should identify at least 
one pole of the T-matrix in the near vicinity, and fits without a nearby pole must be dropped. 

Finally, there is no experimental indication for an S* = — 1, J = 1 s-wave baryon resonance 
below KN threshold, the lowest possible candidates being S(1480) (one star resonance) and 
E(1560) (two star resonance) which are listed as "bumps" in [34]. While spin and parity of both 
states have not been determined yet, recent experiments [42,43] yield controversial results even 
on the existence of E(1480). If, however, these low lying I — 1 resonances should be confirmed 
in the future and have the required quantum numbers, their position is still above the relevant 
energy region considered here. Therefore we drop fits which entail a pronounced 1=1 resonance 
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pole posi1 
first pole 


Gion (MeV) 

second pole 


"WT" 


1420+}^ - 1 (20;«2) 


14401^7 - ^ (76;^}) 


"WTB" 


142311° - ^ (l5;^o) 


1366l^f - 1 [M-f^,) 


"full" 


14181^° - ^ (3l;^3 


1348l^f - z (62;^?,) 



Table 2: Positions of the first and second 7 = pole in the complex W plane. The central 
values correspond to the best fit in each approach, while the error ranges include all fits in the 
1(7 confidence region. 

structure below W = 1.44 GeV caused by an isospin one pole close and immediately connected 
to the real axis. More precisely, fits with 1 = 1 poles located at Im(W^) > —50 MeV, i.e. less 
than twice as far from the real axis as typical A(1405) poles, are not taken into account. 

We observe fits which agree with the DEAR results at a lower overall x^/ d.o.f. than indicated 
in Fig. 4 (but still outside the la confidence region) if the above mentioned criteria are omitted, 
e.g., a x^/d.o.f. = 2.0 value is obtained in the full approach. All of these fits have an isospin 
one pole in common which is very close to the KN thresholds and either a few MeV above (i.e. 
on the physical sheet) or below (i.e. on an unphysical sheet) the real axis. Solutions of this type 
have been reported on in [5,6]. However, such fits clearly violate one of the criteria discussed 
above and are not considered here. 

In the remainder of this section wc will focus on the resonance pole structure of the A(1405), 
i.e., the / = poles that are located on the unphysical sheet which is directly connected to the 
physical real axis between the ttS and thresholds. The nature of the A(1405) has recently 
attracted considerable interest. It has been claimed that instead of the usual appearance of one 
resonance pole the A(1405) results from a pronounced two-pole structure with both poles being 
very close to the physical region [13]. While only the Weinberg- Tomozawa contact interaction 
was taken into account in [13], the inclusion of the next-to-leading order contact terms destroyed 
the pronounced two-pole structure, as one pole was shifted further away from the real axis and 
its contribution to the physical region dissolved in the background [4]. As already pointed out 
in [4] , the position of this pole depends very sensitively on the values of the parameters, whereas 
the other one remains relatively fixed and close to the physical axis. In [5], e.g., which also 
includes the ncxt-to-leading order contact interactions the pole is located at 1321 —i 43.5 MeV, 
i.e., even below the ttS threshold(s) and hence not immediately connected to the physical region, 
whereas in the "WT" approach a pole around a mass of 1390 MeV was found in [13]. 

For most fits in the la confidence interval we observe two isospin zero poles in the region 
ReW = (1250 . . . 1600) MeV, ImW = (-2.5 . . . -250) MeV. In some cases, however, it happens 
that there is only one / = pole in this region. We then extend the pole search beyond the 
chosen limits until a second 7 = pole is found. The observed pole positions are compiled 
in Table 2, where "first pole" refers to the pole which is closer to the real axis at 1.405 GeV, 
i.e., the position of the A(1405) peak. While the variation of the position of this first pole is 
remarkably small in all three approaches, "WT", "WTB", "full", the position of the second 
/ = pole scatters over a wide range in the complex W plane and consequently does not have 
in all fits a significant impact on physical observables. 

We also observe very few fits with a third 7 = pole which appears either at Re > 1.5 GeV 
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and thus well above the KN threshold (s) (recall that the pertinent Riemann sheet is connected 
to the real axis below these thresholds) or deep in the complex W plane (Iml^ < — ISOMeV). 
In these fits, the third pole is thus not expected to have much influence on physical observables. 

From the discussions above it becomes clear that the analytic continuation to the complex 
energy plane and the resulting pole positions depend sensitively on the dynamical input of the 
chiral SU(3) effective Lagrangian. A rigorous extraction of the pole positions, in particular 
the second one, appears therefore very unlikely from the experimental data considered in the 
present investigation. Additional experimental input, however, may help to further constrain 
the position of the second pole, see e.g. Ref. [27]. But in any case, as illustrated above, the 
pole positions can very well serve as an additional constraint to rule out some fits. 

5 Conclusions 

In this work, we have provided an extraction of the strong K~p scattering length derived within 
chiral unitary approaches from fits to available low-energy K~p scattering data. To this end, 
we have utilized three variants of such chiral unitary approaches which differ in the choice of the 
interaction kernel. In the first approach, the interaction kernel is derived from the Wcinberg- 
Tomozawa contact interaction at leading chiral order which is successively supplemented by the 
Born terms and the contact interactions of next-to-leading chiral order in the second and third 
framework, respectively. The usage of three different interaction kernels helps to estimate the 
inherent model dependence of such approaches. 

For all three approaches a least-squares fit to low-energy data in S* = — 1 meson-baryon 
channels is performed. These are in detail K~p scattering, the ttS mass spectrum and the 
precisely measured K~p threshold decay ratios. Fits with the lowest x^/d.o.f. value are found 
in the full approach including the higher order couplings, while the x^/d.o.f. value is largest in 
the Wcinberg-Tomozawa approach. 

Based on a very large variety of different fits to data we can provide an error range for 
the strong K~p scattering length which is related to the strong interaction shift and width 
in kaonic hydrogen. We obtain an energy shift and width in kaonic hydrogen which is in 
agreement with the KEK experiment, but disagrees with DEAR. The present analysis confirms 
the findings of [4] by pointing on questions of consistency of the recent DEAR measurement 
with previous K~p scattering data. The conservative error range for a^-p derived from chiral 
unitary approaches is in clear disagreement with the one deduced from the DEAR experiment. 

Furthermore, we have critically investigated the pole structure of the fits. The first isospin 
zero pole remains relatively fixed in all fits and close to the physical axis, whereas the second 
pole is quite sensitive to the chosen parameters of the approach. In particular, the influence 
of the second pole on physical observables is substantially reduced, if it is further away from 
the real axis, and can even dissolve in the background. Although the pole positions depend 
sensitively on the dynamical input of the chiral SU(3) effective Lagrangian, we have illustrated 
that the general pole structure of a fit can serve as an additional criterion to consider the fit 
as unphysical. In this respect, we look very much forward to the electromagnetic production 
data of the A(1405) from the ELSA accelerator at Bonn which may help to further clarify the 
pole structure of the K~p scattering amplitude below threshold. 
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A Results of the fit to scattering data 

In the appendix, we present the results of the fits to scattering data. The overall agreement 
with the experimental data is very good. In Figs. 5, 6, 7 we plot the elastic and inelastic K~p 
scattering cross sections for the three approaches "WT" , "WTB" , and "full" , respectively, while 
the ttE event distribution is shown in Fig. 8. For the accurately determined threshold branching 
ratios 7, Re, Rn defined in Eq. 6 we obtain: 





"WT" 


'WTB" 


"full" 


exp. [21,22] 
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9 or +0.07 
^•"-•^ -0.06 


9 oe +0.03 
^•"-"J -0.03 


9 oe+O.lO 
^•"-"J -0.09 


2.36 ±0.04 


Rc 


n err +0.001 
U.OOO _o.018 


n 664 +° °22 

U.UU^ -0.024 


U.OOO _o.oi8 


0.664 ±0.011 


R„ 


U.i.Jl _() 


n 1 no +0.009 
U.I.Jo _() 017 


n inn +"-"26 
u.i.JU _o 


0.189 ±0.015 



The central values correspond to the best fit in each approach, while the errors indicate the la 
confidence region. 
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Figure 5: Total cross sections for K^p scattering into various channels calculated in the "WT" 
approach. The best fit is represented by the solid line while the shaded area indicates the la 
confidence region. The data are taken from [15] (empty squares), [16] (empty triangles), [17] 
(filled circles), [18] (filled squares), [19] (filled triangles), [20] (stars). 
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Figure 6: Total cross sections for K~p scattering into various channels calculated in the "WTB" 
approach. The best fit is represented by the solid fine while the shaded area indicates the la 
confidence region. The data are taken from [15] (empty squares), [16] (empty triangles), [17] 
(filled circles), [18] (filled squares), [19] (filled triangles), [20] (stars). 
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Figure 7: Total cross sections for K^p scattering into various channels calculated in the "full" 
approach. The best fit is represented by the solid line while the shaded area indicates the la 
confidence region. The data are taken from [15] (empty squares), [16] (empty triangles), [17] 
(filled circles), [18] (filled squares), [19] (filled triangles), [20] (stars). 
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Figure 8: 7r^S+ event distribution for the three different approaches. The best fit is represented 
by the solid line, while the shaded area indicates the la confidence region. The data are taken 
from [23] and supplemented by statistical errors following [44]. 
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